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Abstract 

The gravitational field of two identical rotating and counter-moving dust 

beams is found in full generality. The solution depends on an almost arbitrary 

function and a free parameter. Some of its properties are studied. Previous 

particular solutions are derived as subcases. 
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I. INTRODUCTION 



The study of the gravitational field of fight beams (null dust) has a long history. In the 
linear approximation to general relativity this was done 71 years ago [1,2]. Later, Bonnor 
found exact solutions which belong to the class of pp-waves (algebraically special solutions 
of Petrov type N and no expansion) [3,4]. He showed that beams shining in the same 
direction do not interact. The gravitational field of two counter-moving light beams is 
very complicated and particular solutions were found only recently [5,6]. One of them was 
generalized to the case of two colliding non-null dusts and named the double-dust solution 
[7]. It is based on a one-parameter relation between two metric components. 

In this paper we find the general double-dust solution and study its properties. It depends 
on one almost arbitrary function and one free parameter. All previous solutions are derived 
as special cases. 

In section 2 the field equations are written and their general solution is given in three 
different ways. The condition of elementary fiatness is imposed, which reduces the number 
of free parameters to one. The restrictions on the arbitrary function, the positivity of 
the energy density and the regular character of the metric are studied. In section 3 some 
particular examples are given such as the solutions of Lanczos [8], Lewis [9] and the two 
Kramer solutions (for null and non-null dust). In section 4 the properties of the 4- velocities of 
the two dust beams are studied. It is proved that the general solution satisfies the dominant 
energy condition. The appearance of closed timelike curves is discussed in the stationary 
version of the Kramer double-dust solution. In section 5 the general interior solution is 
matched smoothly to the vacuum Lewis solution at any distance from the axis. Section 6 
contains a short discussion. 
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II. FIELD EQUATIONS AND THE GENERAL SOLUTION 



We shall work in the stationary formalism when the cylindrically symmetric metric is 
given by [4] 

ds'' = -e^" {dt + Adipf + [e''^ (dx^ + dz^) + W^dip^] (1) 

and there is only radial x-dependence. The static metric may be obtained by complex 
substitution: t —>■ iz, z —>■ it, A —>■ iA. The energy-momentum tensor of the two dust beams 
reads 

Tfj,u = l^ikf^K + 1^-21^^^: (2) 

k^k^ = -1, Z^P = -1 (3) 

where /Xj are the density profiles of the rotating and fiowing beams. The only components 
of their 4-velocities are k^ — V^ — V, k^ — — P — Z. Eq. (3) becomes 

To write the Einstein equations we use the combinations of Ricci tensor components 
utilized in Refs. [10-13]. We shall work, however, with curved mixed components The 
only non-trivial ones are the diagonal components and i?^, i?f . The Einstein equations for 
the corresponding T^^ components give after some rearrangement Eq. (4) and 

k" = -u'^ -B, ^e-^", (5) 

W 

k'^ = u'^ - 5, (6) 

A' = aoW^e-^", (7) 



8vr/.e^(-") = + - + 25, (8) 
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W" 

leTT/iZ^e^^'^-") = (9) 

Here ' is a x-derivative, 2/1 — ni + 1x2 and oq is an integration constant resulting from Eq. 
(7). Units are used with G = c = 1. The equations for the rest components are satisfied 
either identically or when (//i — ^2) Z — 0. Therefore, we accept in the following the relation 
A*i = A*2 = A*- In the case of a single beam, we must put instead Z — 0, which halts its 
motion along the axis. 

There are 6 equations (4-9) for 7 quantities; u, k, W, A, ji, Z, V. Let us take u [x] to be 
an arbitrary function. The equations decouple and starting with Eq. (5) and finishing with 
Eq. (4) the unknowns are found in the above-mentioned order. The function k is obtained 
by integrating Eq. (5). Thus we have 

k' {x) ^ kuu' (x) ^ - (u'^ + b) dx (10) 

and ku may be expressed either as a function of x or m if u' {x) or u' {u) is used. The other 
quantities also have double representations 

W = Wok'e^-fo 'w'^^ = Wokuu'e^-fo^, (11) 

A = aoWo r k'e" 'v-'^'^^dx = aoW^ /" k^e ^0 ^-^""du, (12) 
^0 Jo 



ku 2 



«o -4u n 

~^ {ku ~ 1) 6 " — ky,y,U 



e-2(*^-"), (13) 



7/'p-2(fe-w) 1 

V' = ^-=:^e-^" = l^e-^". (15) 

Eq. (14) shows that ku> 2. The second representation requires the passage to a new radial 
variable u or f — and correspondingly g^x becomes 
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where u' = v! [u). 

One can pass to u in another way, by defining the arbitrary function g [u] = k'. Then 
Eq. (5) becomes a quadratic equation for u' with solutions 

2h (u) = 2u' = -g^ ± ^gl-AB (17) 

and ku — g/h. Eqs. (11)-(15) are easily rewritten in terms of g {u). The restriction ku > 2 
holds iff 
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^e~-F{\), A = e-^\ (18) 



where F is an arbitrary positive function, satisfying F\ > ap/S. A sufficient condition for 
the radical in Eq. (17) to be real is 

F > 2alX, (19) 

which encompasses the previous inequality. 

In the above equations we have set u (0) = 0, so the axis of rotation still coincides 
with the 2;-axis. The indefinite integral in Eq. (11) was replaced by a definite one and the 
integration constant Wq. Another constant k (0) = ko appears when Eq. (10) is further 
integrated. The condition for elementary flatness 

lim ^ (e-^^W^ - e^M^) = 1 (20) 

must be satisfled at the axis. A should vanish there which is accounted for by the lower limit 
of the integral in Eq. (12). W should vanish too. One can see from Eq. (7) that A = o iyV) 
when X ^ and can be neglected in Eq. (20) like in many other cases [11-13]. Therefore 
we obtain the static case condition 

\W\ 

lim I I = 1. (21) 
We inset the expression for W from Eq. (11) and take the limit. Eq. (10) gives 



so that u" (0) < and is finite. Then Eq. (11) shows that W (0) ^ due to u' (0) = 0. Eq. 

(21) yields 

al\Wo\^^e''°. (23) 

Hence, elementary flatness gives a relation between the 3 integration constants without 
affecting the arbitrary functions. The constant ko may be put zero by a coordinate change. 
We won't do this for matching purposes. Thus the general solution depends on one free 
parameter and the function u (x), u' {u) or g {u). The only restriction for u' {u) comes from 
the condition ku > 2. For g (u) this leads to Eq. (19), while u (x) should satisfy 

- (u'"^ + B)dx> 2u'. (24) 

Both sides vanish on the axis, hence, a sufficient condition is 

2u" + u'^ + B<Q or -/V2(ji/2y'>^ (25) 

An important physical requirement is > in some region around the axis, which is 
enough for a interior solution. According to Eq. (13) (0) > always. If k^u < 0, then fi is 
positive everywhere. Otherwise the energy density may become negative at some distance 
from the axis. 

Eq. (16) indicates that the change of coordinates a; — > m is always singular at the origin. 
This may be avoided by introducing a third radial coordinate r, such that u — u (r^) [7] . This 
is equivalent to going back to the first version of the solution by choosing u {x) . However, 
some solutions look simpler when u is the radial coordinate. 

III. SOME PARTICULAR SOLUTIONS 

Let us derive several concrete solutions. For a single beam we must set Z = in order 
to use metric (1) and the field equations above. Then Eq. (14) gives constant m, which is 
set to zero. The solution should be found from Eqs. (4-9). We obtain 
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V^l, W^x, ^aoa;^ (26) 

k = ~ay, Stt/x = ^ale-^'. (27) 

This is the cherished Lanczos solution [8] in comoving coordinates. It represents a cyhnder 
of rigidly rotating dust. In the absence of u (x) it depends just on the parameter qq. 

Another important case is = 0, which, at first sight, should lead to vacuum solutions. 
Two obvious candidates are k = 2u and k = u. Both of them belong to the one-parameter 
series k — {a -\-l)u. This case was solved [7] in the static formulation of the problem. It is 
worth to do this in the stationary frame and see the differences. We must have a > 1 and 
/co = 0. Eqs. (14-15) give directly 

= I , = (28) 

a — 1 a — 1 

Eqs. (11,12) yield the expressions 

W^(a + l)Woff^% (29) 

2(2a + l) ^ ^ ^ ^ 

Eq. (23) gives Wq — ^/a^. It is necessary to express /' through /. Eqs. (5-6) yield 

k' = (31) 

Combining the last three equations gives the result 

4 /.-2_.^\ 

^ 4(2a + l)V-^ ^ )■ ^""^^ 
Positivity requires w < 0, / < 1, < 0. This equation may be integrated 

0.^1^ ^ r (-^) r (1) ^ i\ _ (33, 

2 2a + l 

where ^ = / °'+^ and [p, q) is the incomplete beta function. The formula above can not 
be inverted, but this is not necessary, since we have chosen a new radial variable. Inserting 



Eq. (32) into Eqs. (29), (16) and (13) we find explicit expressions for W^g^u and ji, which 
simphfies to 

2 (o + 1) 

The gravitational field is determined completely. 

In the static formulation 0<a<l,/>l,A;>0. Here we have the opposite. The sign of 
a — 1 is changed correspondingly. The model has two parameters, a (remnant of the arbitrary 
function) and Oq. Kramer puts Oq = 1. Then the central density is bounded, Stt/x (0) < 1/2. 
In fact, it can be as big as we want it, as seen from Eq. (34). In principle, one can choose 
an arbitrary density profile, but the extraction of g (u) out of it is not possible explicitly. 
In Ref. [7] it is asserted that / (x) has no analytical expression. This is true, however, Eq. 
(33) shows that x (/) does have one in special functions. 

Let us look at the case a = 1. The factors a — 1 in Z^, cancel each other, while 
Eq. (4) becomes V'^ = Z^. The two beams turn into null-dust and the metric coincides 
with Eq. (21) (with A = •\/3/4) from Ref. [5] when a — 1 is set effectively to 2, oq = 1 
and ^ = 1 + 3r^/16. This is the proper choice of the final radial coordinate r, instead of 
^ = (1 + r^) ^ in the static case. 

The case a = is a true vacuum case. The static metric becomes simply flat spacetime. 
Here we have 

W^x, f^x'-cl, A^-^ e^-" = l, (35) 

with Co being an integration constant. We should compare this metric to the general Lewis 
solution [9,14,15] 

Wl^x, f^ = lx'---^x''+\ (36) 

^^ = ^^ + 6, e2(A:,-«,)^^i(-'-i), (37) 
nl}L 

which includes 4 constants 6, c, I, n. Obviously, Eq. (35) is a particular case of Eqs. (36,37) 
with n = —1, I = 1, c = Co, b = — Cq ^ and 2co = oq. The solution belongs to the Weyl class 
and has a regular one-parameter metric when > Cq. 



IV. PROPERTIES OF THE GENERAL SOLUTION 



Let us study first the properties of k^^ and P. The 4- velocities are geodesic (no accelera- 
tion) and have zero expansion and shear. The vorticity (twist) vector has been calculated 
with the help of GRTensor for the general cylindrical metric (1) and for the general solution. 
We have — and for 

z _ Qu-2k _ QQ {K - 1) 2fc /ooN 

" 2W^^ " 2{ku-2)'' ' ^^^^ 



w 



. . _±1 y,y, (V _ , ZV - ZVI ^ e. (39) 

For Z'*, w"^ and u?* change sign. These expressions coincide in the Kramer case k — {a-\-l)u 
with Eq. (39) from Ref. [7] after the passage to the static metric is done and one takes into 
account the difference between x and C — P- 

Another useful quantity is F = kH^. It reads for the different metrics 



F = 



1 + 2e2('=-")z2l = ^ = (41) 

J ku-2 a-l ^ ^ 



The proof of the dominant energy condition in Ref. [7] may be hfted to the general solution. 
It depends crucially on the fact that F^ > 1, which is obvious from the above formula. 

It is well known that closed timelike curves (CTC) exist when g^^ < 0. This condition 
is rather intractable further for the general solution. In the Kramer case we have 

9.. - ^i^/-^ [(2a^ + 2a + l) - aV"^ - (a + 1)^ f] ■ (42) 
(2a+l) '-^ ^ 

The sign is determined basically by the competition of the first and the second terms in the 
square brackets. It becomes negative when 

e>l + -^. (43) 
(a + 1) 



This happens always for big enough r. In the static metric formulation g^,^^ is strictly positive 
and there are no CTC as seen from Ref. [7], Eq. (36). In this case a matching can be done 
at some tq to the Levi-Civita static metric and a realistic global solution constructed. In 
the stationary case the matching should be done to the Lewis solution (36,37) whose Weyl 
class is essentially static, but the Lewis class contains CTC. 



V. THE GLOBAL SOLUTION 

We shall show that the general interior solution can be matched smoothly to the Lewis 
solution (36,37) for any value Xq of the radial coordinate. For this purpose we perform 
scaling on two of the coordinates in the exterior, i = rt and (p = fl(p. Instead of scaling z we 
introduce again ko in the interior. In Ref. [7] the transformation P — h (p) was also applied, 
where h is an arbitrary function of the exterior radial coordinate p. The general solution, 
however, aheady includes an arbitrary (modulo some positivity requirements) function u {x) 
and it is redundant to introduce a second one. The effect of such transformation then is that 
u{xq) and u' {xq) become free parameters, while quantities given by integrals, like /c(xo), 
are fixed by the interior solution. The matching procedure is the following: we choose an 
interior solution which fixes ao, k (xq), k' (xq), W (xq) , A (xq) and A' [xq) up to Wq. Next 
we demand that the metric should be continuous together with its first derivative at the 
junction. This gives several algebraic equations for the free constants u (xq), u' (xq), fco; 'T, 
n, n, I, c and h. Finally, we fix them by solving the equations, the only free parameter 
remaining being xq. We consider the first representation of the general solution and all 
functions in the following are taken at the point Xq. 

The continuity of the metric yields four conditions 

fir' = /^ flALTn = fA, (44) 
(flAl + fE'Wl) = fA' + r'W\ (45) 
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The continuity of the metric derivatives supphes another four. After some rearrangement, 
the total system of 8 equations becomes 



g2feo^4("'-l)g2„-2fe^ (47) 



- 1 = 4x0 {k' - u') , (48) 



= ^A'„ (49) 



~ = 1~' rn = —, 50) 



w 

= 1, (51) 



flA'^ = aoxof-', (52) 
u'. (53) 



Eq. (47) determines kg. Then Eq. (23) determines Wo and, hence, W. Eq. (48) determines 
n, while Eqs. (37,49) fix b. Eq. (50) yields expressions for Q and r 



T = 



Eqs. (51-53) form a system for /, u' and c/l 



u- = (55) 



c nao^o " 



/ P{l + n-2u'xo) 



(56) 



l-n-2W=^^^. (57) 



Replacing Eq. (56) into Eq. (57) and using Eq. (48) we obtain 
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(58) 



The combination of Eqs. (55,58) gives 



/ = 1, 



u 
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Ak' 



+ Xoal 
4xo 



(59) 



Inserting these formulas in the previous equations, everything is determined in terms of the 
interior solution and Xq. For example, we have for the parameter n, which measures the line 
mass density 



There are several interesting features of the matching. The gravitational potential / does 
not depend on any parameters at the junction and has the same value as at the axis. Only 
the ratio of I and c is determined. The quantity is not necessarily positive. When < 
we enter the Lewis class of the Lewis solution, which possesses CTC [14,15]. 



We have obtained in this paper the general global stationary cylindrically symmetric 
solution for the gravitational field of two identical, rotating and counter-moving dust beams. 
Three representations of the interior solution have been given. They depend on the free 
parameter and the arbitrary function u {x), u' {u) or g {u), which satisfies the condition 
(25), A;„ > 2 or (19) respectively. A particular solution with qq — 1 and depending on the 
arbitrary parameter a has been found by Kramer [7]. Many of its nice properties are shared 
also by the general solution. It satisfies the dominant energy condition. The energy density 
of the beams is non- negative. The axis is regular and elementary fiat. The solution is 
necessarily non-diagonal. We have studied its stationary alternative. Rotation compensates 
gravitational attraction and prevents collapse and appearance of singularities. 

It was shown that two non-zero components of the dust four- velocities are enough for a 
solution with arbitrary density profile. In the case of colliding null-dust three such compo- 
nents are necessary [13]. 




(60) 



VI. DISCUSSION 
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The interior solution can be matched at any distance xq to an exterior vacuum stationary 
solution, the Lewis solution [9,14,15]. Thus a global solution is formed. An important 
property is the traditional appearance of CTC in rotating cylindrically symmetric solutions. 
This happens both in the interior when xq is big enough and in the exterior, when the Lewis 
class is induced. The Weyl class is locally equivalent to the Levi-Civita solution [14]. It is 
causal and serves as an exterior for the general static solution. 
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